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Abstract. We recently proposed in a Letter [Physical Review Letters 108 255303] a novel 
scheme to detect topological edge states in an optical lattice, based on a generalization 
of Bragg spectroscopy. The scope of the present article is to provide a more detailed and 
pedagogical description of the system - the Hofstadter optical lattice - and probing method. 
We first show the existence of topological edge states, in an ultra-cold gas trapped in a 2D 
optical lattice and subjected to a synthetic magnetic field. The remarkable robustness of 
the edge states is verified for a variety of external confining potentials. Then, we describe 
a specific laser probe, made from two lasers in Laguerre- Gaussian modes, which captures 
unambiguous signatures of these edge states. In particular, the resulting Bragg spectra 
provide the dispersion relation of the edge states, establishing their chiral nature. In order 
to make the Bragg signal experimentally detectable, we introduce a "shelving method", 
which simultaneously transfers angular momentum and changes the internal atomic state. 
This scheme allows to directly visualize the selected edge states on a dark background, 
offering an instructive view on topological insulating phases, not accessible in solid-state 
experiments. 

1 Introduction 

Ultra-cold atoms in highly controllable coupling fields constitute a novel experimental tool for studying 
the rich many-body physics arising in two dimensions [1,2,3,4]. Motivated by the possibility of reaching 
interesting quantum phases, synthetic magnetic fields [5] and spin-orbit couplings [6,7,8] have been 
realized experimentally for neutral atoms. Today, the engineering of these synthetic gauge potentials 
opens an important path for the exploration of topological phases, such as quantum Hall (QH) states, 
topological insulators and superconductors, in the clean and versatile environment offered by cold-atom 
setups [2,3]. 

For the last decades, these topological phases have gained the interest of the scientific community 
for their unique properties, such as quantized conductivities, dissipationless transport and edge-states 
physics [9,10]. These impressively robust phenomena rely on an important concept, the so-called bulk- 
edge correspondence [11,12]. Topological phases of matter are characterized by robust, integer- valued 
topological invariants related to the bulk structure of the material. The bulk-edge correspondence stip- 
ulates that well-defined edge excitations localized near the boundaries of the system are associated to 
these topological invariants. Such edge excitations are of tremendous practical importance, as they usu- 
ally carry some form of current protected against perturbations as long as the topological structure is 
preserved. As such, they are at the origin of the dissipationless transport observed for these phases. 
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For example, in the QH effect taking place in 2D electronic systems, the topologically invariant Chern 
number [13,14] guarantees the presence of current-carrying edge states, and imposes their chirality [11]. 

Cold-atom realizations of topological phases therefore constitutes a complementary, but also intrin- 
sically appealing, playground to further deepen our understanding of these topological properties. How- 
ever, the detectability of topological phases remains a fundamental issue in the cold-atom framework 
[15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35], where transport measurements consti- 
tute a possible, but very challenging task [36]. In this sense, alternative signatures of topological phases, 
together with novel experimental probes, have to be considered in this new context. Following this 
strategy, several schemes have been described to directly measure topological invariants, based on spin- 
resolved time-of-flight [27,33] and density measurements [19,28]. Alternatively, Bloch oscillations could 
also be performed to evaluate the Berry's curvature in 2D atomic systems [30], which could then provide 
an estimation of the Chern number when integrated over the Brillouin zone. 

Inspired by the bulk-edge correspondence, it has also been suggested that topological edge states 
could be directly probed [21,22,37,29]. For example, in the context of cold-atom QH insulators, a sat- 
isfactory signature of the non-trivial topological order would be obtained by probing the dispersion 
relation of QH edge states, thus demonstrating their chiral nature. 

It is the aim of the present work to describe in detail such a realistic probe. We choose to analyze this 
detection scheme for an optical-lattice setup reproducing the Hofstadter model [38], which is one of the 
simplest tight-binding lattice model exhibiting non-trivial Chern numbers [13,14] and topological edge 
states [11]. The experimental realization of this model, using cold atoms in optical lattices, is currently 
in development in several laboratories [39,40,41], based on the proposals [42,43]. We believe that our 
detection scheme could easily be extended to any ultracold-atom setup emulating 2D topological phases. 



2 The Hofstadter optical lattice and topological edge states 

We start with a two-dimensional fermionic gas confined in a square optical lattice and subjected to a 
uniform synthetic magnetic field B = Biz [42,43]. The Hamiltonian is taken to be 

~ ~ J ^m+l,n^^)^ ^ ^ ^m,n+l^^)^ ^ h..(Z. + ^onf(^) ^\n,n^rn,ni (1) 

m,n m,n 

where Cm,n is the annihilation operator defined at lattice site (m, n) G 1? and where J is the tunneling 
amplitude. The site indices m, n = 1, . . . are related to the spatial coordinates through 

(x, y) = a{m - N/2, n - N/2), 

such that the center of the system is defined at xq = (0,0). In the following, the lattice parameter a 
defines our length unit. The external and cylindrically symmetric confining potential is written as 

Konf(r) = J(r/redge)^ (2) 

where 7 = 2 (7 = 4) corresponds to the standard harmonic (quartic) trap used in cold-atom experi- 
ments. The expression used for the confinement (2), including the tunneling amplitude J, is chosen such 
that the edge of the atomic cloud - the Fermi radius - is given by ~ Tedge, for the specific config- 
uration considered in this work (i.e. ^ = 1/3 and Ep ^ —1.5 J, cf. Section 2.2). In the absence of the 
confinement Konf = 0, Eq. (1) describes the Hofstadter lattice model, namely, a gas of non-interacting 
fermions in the tight-binding regime, subjected to a vector potential A = (0, Bx, 0) corresponding to ^ 
magnetic flux quanta per unit cell [38]. Methods to implement this optical-lattice setup, illustrated in 
Fig. 1 (a), have been proposed in Refs. [42,15,43,40], and some important experimental steps towards 
this goal have already been successfully achieved [39,40,41]. 
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Fig. 1. (Color online) (a) The Hofstadter Optical lattice: cold atomic fermions are trapped in a square optical 
lattice and subjected to a synthetic uniform magnetic flux (b) The corresponding bulk energy spectrum 
E = E{k) for ^ = 1/3, and in the absence of external confining potential. The Chern numbers Ni^2 associated 
with the two lowest bulk bands, and computed through Eq. (3), are shown. The Hall conductivity gh^ computed 
from Eq. (4), is indicated for Fermi energies located inside the two bulk gaps. 
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Fig. 2. (Color online) The cylinder analysis, (a) the Hofstadter optical lattice with ^ = 1/3 defined on an 
abstract cylinder C, in the absence of external confining potential. The two opposite edges of the cylinder are 
designated by OCt and OCb, and an arbitrary unit plaquette is emphasized, (b) The corresponding energy 
spectrum E = E(ky), as a function of the quasi- momentum. The three bulk bands are designated by purple 
shades (see also the bulk spectrum in Fig. 1(b)). The dispersion branches located inside the bulk gaps correspond 
to chiral edge states, which are spatially localized at the two boundaries OCt^b of the cylinder. The edge states 
associated with different gaps have opposite velocities Ve = E/hky. 



2.1 The topological edge states and the bulk-edge correspondence 



The transport properties and topological phases exhibited by the system can be deduced by solving 
the single-particle Schrodinger equation, which is performed through a direct diagonalization of the 
Hamiltonian (1). Setting Konf = and considering periodic boundary conditions along the x and y 
directions - namely, solving the system on a two-dimensional torus T^- one obtains the energy band 
structure E = E{k) depicted in Fig. 1 (b), where k = ikx^ky) is the quasi-momentum. For ^ = p/q^ 
where p^q e the spectrum splits into q subbands [38], separated by bulk gaps. In the following, we 
set <P = 1/3, in which case the spectrum depicts two large bulk gaps of the order A ^ J. When the 
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Fermi energy is set within a bulk gap, e.g. Ep ^ ±1.5J, the interior of the system has the properties of 
an insulator. 

Because of the chosen periodic boundary conditions, the toroidal geometry cannot account for the 
edges present in real systems. In most physical systems, edge effects are normally neglected for suffi- 
ciently large sample size. However, topological phases constitute a counter-example, where the behavior 
at the edges represents an essential component of the physics. To see this, let us consider the same model 
on a cylindrical geometry, where periodic boundary conditions still hold in the y direction, but where the 
system has a finite length in the x direction. In this geometry, which is still abstract from the cold-atom 
point of view, the lattice system now features two edges, as represented in Fig. 2 (a). The corresponding 
energy spectrum E{ky)^ shown in Fig. 2 (b), can be partitioned in terms of hulk states and topological 
edge states. Indeed, one finds two dispersion branches within each bulk gap (cf. light blue and red curves 
in Fig. 2 (b)), describing states that are spatially localized at the two edges of the cylinder [11]. Note 
that the dispersion relation is approximately linear in the middle of the gap, E/h ^ v^ky where is 
the group velocity, a feature generically found for such models [10]. In particular, we find that the edge 
states located in different gaps have opposite chirality sign('Ue). Therefore, when the Fermi energy is set 
within a given bulk gap, low-energy excitations propagate along the edge of the system with a specific 
chirality. In the solid-state framework, these chiral states are responsible for the quantum Hall effect 
[9,10,11]: when E^ ^ ±1.5 J, this propagating edge structure leads to the quantization of the transverse 
Hall conductivity [11], cfh = ±1 (in units of the conductivity quantum), as the two bulk gaps host a 
single edge-state branch per edge, with opposite chirality sign('Ue). 

In fact, these chiral edge states are topological^ in the sense that they are directly related to topological 
invariants associated with the bulk gaps. This important fact guarantees their robustness against small 
external perturbations: the edge states survive as long as the non-trivial bulk gap in which they reside 
remains open. The fundamental concept which relates the edge states to topological invariants is the 
so-called hulk-edge correspondence [11,12], which is briefly described below. First, let us introduce the 
Chern numbers Nj^^ which are topological indices defined for each bulk band Ejy{k), labeled by the index 
z/, through the Thouless-Kohmoto-Nightingale-Nijs expression (TKNN) [13], 



This formula corresponds to the integral over the first Brillouin zone (T^) of the Berry's curvature 
associated with the eigenstate \ujy{k)) belonging to the band Eiy{k). Here, k = (kx^ky) G is the 
quasi-momentum. When the Fermi energy is exactly located in a bulk gap, the Hall conductivity is 
directly related to the Chern numbers. 



which can be derived from the Kubo formalism [13,14]. Here the conductivity is expressed in units of 
the conductivity quantum. For the Hofstadter lattice (1) with ^ = 1/3, that we consider in this work, 
the bulk energy spectrum splits into three energy bands, which have the associated Chern numbers 
A^i = — 1, A/'2 = 2 and = —1 (cf. Fig. 1 (b)). Therefore, when the Fermi energy lies in the first 
(second) bulk gap, the Hall conductivity corresponds to an = —'^ {cfh = +1), as illustrated in Fig. 1 
(b). In this sense, the Hall conductivity gh = ±1 is a non-trivial topological index characterizing the 
two bulk gaps ^. 

In such a non-trivial topological band configuration, the bulk-edge correspondence dictates the fol- 
lowing result [11,12]: if we solve the model (1) on an open geometry, such as the cylinder considered 
above, gapless edge states will necessarily appear in the bulk gaps, because the latter are associated 
with non-zero topological invariants. Moreover the number of edge-state branches inside a bulk gap is 

^ The Hall conductivity associated with a bulk gap can only change its value through a topological phase 
transition, that is, a gap-closing process. 




(3) 




(4) 
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given by the modulus of the Hall conductivity in (4), and their chirality by sign{aH)' This fact is easily 
verified by comparing the numerical results shown in Figs. 1 (b) and 2 (b). The presence of a single 
edge excitation per physical edge in the lowest bulk gap of Fig. 2 (b) is in agreement with the fact 
that I cTj:/ (lowest gap)| = |A^i| = 1: this is precisely the bulk-edge correspondence applied to the present 
context. 



2.2 The Hofstadter optical lattice in an external confining potential 

Having analyzed the edge-state structure, based on the cylinder analysis presented above, we now come 
back to the actual optical-lattice setup, whose finite size is determined by the confining potential Konfl^)- 
For an infinitely abrupt potential (7 = oo), 

Vconf{r < Tedge) = 0, KonfC^ > ^edge) = OO, (5) 

the lattice has the disk geometry represented in Fig. 3 (a). The corresponding spectrum and eigenstates 
are shown in Fig. 3 (b) and Fig. 4 (a). From the latter result, we find that the clear partition of the spec- 
trum in terms of bulk states and topological edge states still holds in the experimental planar geometry: 
similarly to the cylindrical case, we obtain bulk states within the bulk bands and edge states (illustrated 
in Fig. 3(d)) within the bulk gaps. This result is in agreement with the bulk-edge correspondence, which 
requires that the chiral edge states lying inside the bulk gaps do not depend on the particular geometry 
of the lattice, as they are dictated by the Chern numbers (3) associated with the bulk. Therefore, when 
considering the realistic circular geometry produced by the confining potential Konf(^) in Eq. (2), one 
obtains the same edge-state structure propagating along the circular edge r = redge as the one obtained 
from the abstract cylinder discussed above: the number of edge- state branches (per physical edge) and 
the chirality deduced from them are identical, as these properties do not depend on the chosen geome- 
try. In other words, the only difference between the cylindrical and the circular Hofstadter model is the 
number of physical edges (i.e. two edges for the abstract cylinder, and only one edge for the realistic 
circular geometry). However, let us comment on the fact that real boundaries, with finite 7 7^ oo, do 
slightly affect the dispersion relations - and thus the angular velocity - of the edge states (cf. paragraph 
below). 

The bulk-edge correspondence indicates that the lowest bulk gap in Fig. 3 (b) hosts a single edge- 
state branch with a negative angular velocity, since the corresponding Hall conductivity an = — 1 is 
solely governed by the topological expression (4). Furthermore, the edge-state branch present in the 
second bulk gap corresponds to the opposite chirality, since a// = +1 when the Fermi energy is in the 
highest gap. These results have been verified by directly computing the angular velocity of the edge 
states, using the expression 

Oe = {i/ h) ^ V^* (m, n) [fiQ , 0] i/je {m,n), 

m,n 

= ^ ^(6>(m + 1, n) — 6>(m, n))ilj2{m^ n)'il)e(jn + 1, n) 

m,n 

+ (^(m — 1, n) — ^(m, n))tlj*{m^ n)tlJe{Tn — 1, n) 

+ e^^^^^ ((9(m, n + 1) - 6>(m, n)) V^*(m, n)V^e(m, n + 1) 

+ e-^^^^^((9(m,n- 1) -6>(m,n))V^*(m,n)V^e(m,n- 1), (6) 

where ^(m, n) = atan2((n — A^/2), (m — 7V/2)) and m, n = 1, . . . , A^. The numerical results Oq ^ ±0.07 J/h 
are shown in Fig. 3 (b), together with a sketch of the corresponding dispersion relation in Fig. 3 (c). 

The partition of the energy spectrum in terms of bulk and edge states remains valid for finite 
confinements [31]. This fact is demonstrated in Fig. 4(a)-(d) for different values of the parameter 7 = 
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(a) (b) (c) 




Fig. 3. (Color online) (a) The Hofstadter optical lattice in the cylindrically-symmetric confining potential 
V^onf(^), with 7 = (X) and redge = 13a. A unit plaquette is highlighted on the figure, (b) The corresponding 
discrete energy spectrum where the index a classifies the eigenvalues in ascending order. The energy ranges 
corresponding to the bulk states are designated by purple shades (cf. Fig. 1). The angular velocity {0) associated 
with the edge states at e = ±1.5 J, computed through Eq. (6), is shown in the gaps, (c) Sketch of the dispersion 
relations corresponding to the edge states around e ^ ±1.5 J. At zero temperature, all the states are occupied 
below the Fermi energy, (d) The amplitude \il:e(Tn^n)\^ corresponding to the edge state at Ce ~ — 1.5J = £^f, 
highlighted by a purple dot in (b). (e) The probe shape //.(m, n) ^/^(r) used to detect the edge state in (d). 

2,4, 10, defined in Eq. (2). In this figure, we observe how the edge-state structure smoothly follows the 
Fermi radius i^p, i-e., the edge of the atomic cloud imposed by the external confinement. In particular, 
for Ey ~ —1.5 J, we find that the edge states remain close to the parameter redge ^ defined in 
Eq. (2). In the following, we will consider a Fermi energy = —1.5 J located within the first bulk 
gap, so that the fermionic gas forms a QH insulator, with central density n = l/3a^, and such that the 
edge states are located close to the radius redge- Let us finally emphasize that the edge states velocity 
highly depends on the boundary produced by the confinement: Oq significantly decreases as the potential 
Konf(^) is smoothened, e.g.^ Oq ~ — 0.02 J/h for 7 = 10 and Oq ~ —O.OlJ/h for 7 = 2. 



3 Angular Momentum Spectroscopy : 

In the previous Section, we showed that the two bulk gaps at e ~ ±1.5 J host topological edge states with 
opposite chirality. The core of our proposal is to design an experimental probe, yielding a clear signature 
from these topological states, exploiting their specific chirality to distinguish them from the bulk. This 
probe is inspired by Bragg spectroscopy [22,21], a form of momentum-sensitive light scattering which 
is routinely performed to access the linear momentum distribution of cold atomic gases [44,45]. First of 
all, we note that for the present problem, in which chiral edge states propagate along the circular edge 
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(a) 7 = oc |^,(x,, = o)p (b) 7 = 10 mx^y = ^)\' 




Fig. 4. (Color online) Confinement effect: The amplitudes |'0e(a:,2/ = 0)|^ associated with the single-particle 
states iV^e) of the confined Hofstadter optical lattice (cf. Fig. 3 (a)-(b)), as a function of their energy e and spatial 
coordinate x. Here the y coordinate is chosen to be at the center of the trap, y — 0. The confining potential 
is taken to be yconf(^) = J{r/redge)^ , with redge = 13a, and (a) 7 = oo, (b) 7 = 10, (c) 7 = 4 and (d) 7 = 2. 
The first bulk gap around Ef — —1.5 J, indicated on the figures, correspond to the energy range where chiral 
edge states are found. The radius r = redge, at which the Bragg probe is focused to scatter the edge states, is 
indicated by white vertical lines. Note that the Fermi radius Ry ~ Tedge, for Ey ~ —1.5 J. 



of a 2D disk (cf. Fig. 3 (a)), it is more convenient to probe the angular momentum distribution in the 
vicinity of this edge^. Therefore, we propose 



^ A standard Bragg spectroscopy, measuring the linear momentum distribution, could be used to probe the 
system as well. However, to achieve good overlaps with the edge states, and to avoid spurious signal from the 
bulk states, the Bragg lasers would have to be focused on a small region near the edge of the cloud, where one 
of the velocity's component, say Vx^ remains approximately constant. Using simple dimensional arguments, one 
expects an extremely small excitation rate in this case. In contrast, our scheme exploiting angular momentum 
allows one to focus the probe lasers on the entire radius redge, maximizing the overlap with the edge states 
spatial distribution and thus the excitation rate. 
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— to use a spatial mode carrying angular momentum, in order to probe the angular momentum distri- 
bution; 

— to shape the probing lasers to maximize (minimize) the probability to excite edge (bulk) states. 

We consider two lasers in high-order Laguerre- Gauss modes, denoted 1, 2, with optical angular momenta 
^1,2, which correspond to the electric fields 

^i,2(r) = y^2fi,,iAr)eM-iK20 - ioJi,2t), fi{r) ^ (r/ro)l^le-^'/2^o , (7) 

where (r, 0) are polar coordinates. The beams are assumed to be set off-resonance from a neighboring 
atomic transition, so that spontaneous emission can be neglected. This leads to a scattering Hamiltonian 

^BraggW = hf2^{ll^e-'--'+Qe'--')cicp, (8) 

11^ = 1 1 dxraix)M^)fLir)e"^', (9) 

where the index q = l2 — h represents the amount of angular momentum transferred by the probe (in 
units of h), huL = h{uji — UJ2) is the energy transfer, and O is the Rabi frequency characterizing the 
strength of the atom-light coupling. The probe profile is 

/L(r) = (r/ro)^e-'-'/'-o/AAi, (10) 

with L = |/i| + 1^2 1, and is illustrated in Fig. 1 (d). The operator cj^ creates a particle in the eigenstate 
li/ja) of the unperturbed single-particle Hamiltonian, i.e. ho\ilJa) =^a\i^a)' Let us rewrite the scattering 
Hamiltonian as 

^Bragg(t) = hf2{W,e-'"-' + W^_,e--*), W, = Y.Ilpcicp. (11) 

q;/3 

Solving the time-dependent problem Ho-\-HBragg{t) to first order, we write the many-body wave function 
as 

\m) = Mm + Yl ^mW^-^^^'/^i^m) ^ Mm + E hi{t)e-^^'^'\ki), (12) 

M {k,l) 

where |0) = ni/<£;F ^t\^) denotes the groundstate at zero temperature, and 

\kl) = \1... 1^1... 1 I 0...0^0...0), (13) 

where k > E^^ I < Ey and uoki = (e/c — ei)/h > 0. Here we have restricted the full Hilbert space to the 
subspace spanned by the ground state and the excited states that are coupled to it to first order in the 
perturbation (11). 

Setting the initial condition (6o(0) = ^^bki{0) = 0), one finds 

b,i{t) = -if2(^4Vfe-i('^L)e^^-* + (J,l)*5+(c.i)e*^^<*), (14) 

where S^i{ujl) = sin(Z\^^t)/Z\^^, zA^^ = {ujki ±a;L)/2. The number of scattered atoms, or excitation 
fraction, is then given by 

N{q,u;L) = Yl I^'^'WI' = E + Wfc)*5+(c^L)e*<'?. (15) 

k,l k,l 
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In the long-time limit, and neglecting the anti-resonnant term (oc e*^fc^^), this yields the standard Fermi 
golden rule 

N{q,ojL) =27r f2H ^ \qfS<^'\uki-WL), (16) 

k>EF,l<EF 

t^oo 

where 6^^\uj) = (2/7rt)sin {ujt/2)/uj'^ > S{uj). The expression (16) emphasizes the explicit relation 

between the excitation fraction N{q^ujL) and the rates defined in Eq. (9). When the Fermi energy 
Ep is set within a bulk gap, and for small frequency and intensities ujl-,^ ^ J /^^ the excitation fraction 
N{q^ujL) probes the dispersion relation ee = ee(M) associated with the gapless edge states |?/^e) that lie 
within this gap, where M is a quantum number analogous to angular momentum (see Fig. 5 (a)). For an 
optimized probe shape /l(^), which we obtain by setting ro ~ ^edge/ this can be deduced from 
the behavior of the overlap integrals /^^ defined in Eq. (9). They are represented in the uo^i — Q plane 
in Fig. 5(b) for 7 = oo. At low frequencies Uki <C J/h^ we find a continuous alignment of resonance 
peaks ujj^i^ ~ OeQ. This reflects the linear dispersion relation e/h ^ OqM in the vicinity of the Fermi 
energy, and provides the angular velocity Oq ^ — 0.07 J/h and the chirality (i.e. sign(^e)) characterizing 
the edge states in the lowest bulk gap. We find that this result is in perfect agreement with the direct 
evaluation of the angular velocity, obtained through Eq.(6). As already mentioned, we stress that the 
edge states velocity significantly decreases as the potential Fconf(^) is smoothened. The absence of sub- 
stantial response for > in Fig. 2 (b) clearly proves that our setup is effectively sensitive to the edge 
state chirality. Naturally, the signal obtained by setting the Fermi energy in the second bulk gap, or by 
reversing the sign of the magnetic flux ^ — would probe the opposite chirality. 

At finite times, it is preferable to evaluate the excitation fraction through a numerical evaluation of 
the Schrodinger equation, 

.^db^ = hnJ2Wki;uUt)buUt)e'^^'"-''-^'^\ (17) 

n,m 

where Wki;nm{t) = {kl\W {t)\nm) and W{t) = #<^e"^^^^ + W'_ge^^^^ The many-body wavefunction is 
still restricted to the first-order subspace but off-resonant terms and deviations from the long-time limit 
are included. For the reasonable finite times and small Rabi frequencies Q <^ J/h used in our calcula- 
tions, we find that the excitation fraction N{q^ujL) obtained from a numerical evaluation of Eq. (17) is 
in perfect agreement with Eq. (15). A typical result is presented in Fig. 6, for q = —4 and = —1.5 J, 
emphasizing the three distinct regimes of light scattering: "edge-edge" , "bulk-edge" and "bulk-bulk" . 
The "edge-edge" regime corresponds to transitions solely performed between the edge states close to 
Ep: A sharp resonance peak is visible at uj^£^ ~ Oeq ~ 0.3 J/h for q = —4, and stems from four transi- 
tions between edge states, as sketched in Fig. 5(a). Then, at higher frequencies, ujl ~ J/K small peaks 
witness allowed transitions between the lowest bulk band and the edge states located above E^. Finally, 
for ujl ^ 2J/h^ many transitions between the two neighboring bulk bands lead to a wide and fiat signal. 
As shown in Fig. 7 (a), this bulk-bulk response is significant for both q = ±4, as a consequence of 
the large density of excited states in this frequency range. In the following, we consider the quantity 
N{q,ujL) — N{—q^ujL)^ which is zero for a system with time-reversal symmetry (cf. Fig. 7(b)). We have 
repeated the calculations for several potential shapes, finding no qualitative change (cf. Fig. 3(b)). Al- 
though it is advantageous to use a steep confining potential, where the edge states are exactly localized 
at r = Tedge, the signal from the edge states is robust even in a harmonic trap (7 = 2). We stress that a 
well focused probe allows to significantly reduce any signal from the bulk. Finally, we note that excita- 
tion times of several h/J^ which seem experimentally realistic, are long enough to resolve the edge-edge 
resonance but still too short to neglect the broadening due to the finite pulse time. 

As can be deduced from Fig. 5(a), the number of allowed transitions \ipk) scales with the 

probe parameter q in the "edge-edge" regime. This is a consequence of the linear dispersion relation in 
the vicinity of Ep. Thus, one observes an increase of the peaks for increasing values of \q\ (cf. Fig. 7(c)). 
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0.1 0.2 0.3 . 0.1 0.2 0.3 

Fig. 5. (Color online) (a) Sketch of the single-particle energy spectrum and the transitions {ipi) \ipk) 
probed by N^q^uol)^ for uol ~ (jOki ^ J/h. (b) The amplitude l/^J^, as a function of the probe parameter q 
and excitation frequency ujki^ for # = 1/3, Ey — —1.5 J, L = 13 and ro = 5.1a. (c) The amplitude as a 

function of the probe parameter q and excitation frequency ujku for Ey — —1.5 J and Ey — +1.5 J. The confining 
potential is infinite (7 = 00) and redge = 13a in all figures. The angular velocity of the edge states present in 
the first [resp. second] bulk gap is 0^ ^ (^IT/q ~ —0.07 J/ h [resp. Oe ^ -\-0.07J/h]. Thus, the two bulk gaps are 
associated with opposite chiralities, in agreement with the result obtained from Eq. (6). 



We stress that this progression only occurs in the "edge-edge" regime, namely when q is chosen such that 
huj^£^ ~ hOeq is smaller than the energy difference between and the closest bulk band. In the case 
illustrated in Fig. 7(c), the "edge-edge" regime is delimited by \qe-e\ ^ 7. Beyond Ig'e-el, the resonance 
peak enters the "bulk-edge" regime. In this case, the excitation fraction N{q^ujL) broadens, N{—q^ujL) 
is no longer negligible, and the linear dispersion relation is no longer probed. We thus conclude that 
a moderate value (here \q\ ~ 4) is preferable to keep a narrow peak, well separated from the broader 
"edge-bulk" signal. 



4 Isolating and imaging the edge states: The Shelving method 

In the previous Section, we showed that the Bragg spectra allow to reconstruct the dispersion relation 
of the edge states. Thus, from this, we could naively claim that the Bragg probe (11) offers a clear 
signature of the edge states, suitable in a cold-atom environment. However, from an experimental point 
of view, the Bragg scheme described above presents a major drawback. Indeed, the number of scattered 
particles N{q^ujL) can only be measured through the modification of the spatial or momentum densities 
after a Bragg pulse, and unfortunately, these signatures only constitute small perturbations on top of the 
strong "background" of unperturbed atoms. For a circular system with Fermi radius Rf^ one can expect 
about A/edge A/hOe ^ RpA/hlvel edge states, with A the bulk energy gap and the group velocity 
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(a) (1) edge-edge (2) bulk-edge (3) bulk-bulk 




= -1.5J 



2 4 



Fig. 6. (Color online) (a) Excitation fraction N{q^ujL) versus probe frequency, for an angular momentum 
transfer q — —4. Here, Q — 0.05 J/h, t — 20h/ J ^ L — 13, ro — 5.1a, redge — 13a, # — 1/3, Ey — —1.5 J and 7 — 00. 
The panels (l)-(3) sketch the three regimes of light scattering, shown in the main figure, for q= —4: (1) Processes 
between edge states located in the gap around Ey — —1.5 J. (2) Processes between the lowest bulk band and the 
edge states located above the Fermi energy. (3). Transitions between the two lowest bulk energy bands. 




^L{J/h) UL{J/h) UJLiJ/h) 



Fig. 7. (a) Excitation fraction N{q,ujL) versus probe frequency, for an angular momentum transfer ^ = ±4. 
(b) N(q,ujL) — N(—q,ujL) for several shapes of the trapping potential V^onf(^) = ^(^/^edge)^ and q= —4. (c) 
Increasing then broadening of N {q^ uj l) — N {—q^ uj l) ^ for increasing |^|. In all the figures i7 = 0.05J//i, t = 20h/J, 
L=13, ro = 5.1a, redge = 13a, ^=1/3, ^f = — 1-5J and 7 = cxo (except in (b)). 

(cf. above). Using the parameters from Fig. 6(a), one finds A^edge ~ 20 while the total number of atoms 
in the calculation is N ^ 200. In two dimensions, we note that the Fermi radius Rp ~ a/^V, so that for 
more realistic particle numbers {N ^ 10^), one can expect about ~ 150 edge states. Detecting these 
few edge states on top of the signal coming from the ^ 10^ unperturbed ones constitutes a significant 
experimental challenge with present-day technology. One possibility to avoid this difficulty is to use an 
alternate detection scheme, where the probe also changes the atomic internal state. The probe signal 
can then be measured against a dark background (without unperturbed atoms), which allows powerful 
imaging methods to be used {e.g. large aperture microscopy, as recently demonstrated for quantum gases 
in optical lattices [46,47]) . 



4.1 Shelving detection method using state-changing transitions: example for ^'^^Yb atoms 

We present in Fig. 8 a possible detection scheme suitable for two-electron atoms with ultra-narrow op- 
tical transitions, inspired by the electron shelving method (in the following, we consider ^^^Yb atoms to 
give a specific example). The ground g and metastable excited e states have zero electronic angular mo- 
mentum but nuclear spin 7 = 1/2. We denote the Zeeman manifolds {g^^g^} and {e^, e^} in the ground 
^5*0 and ^Pq excited states, respectively. The states g^ and are initially populated, as laser coupling 
between these two states is used to generate the artificial gauge field [42,43] leading to Eq.(l). The strat- 
egy, which is at the core of our detection scheme, is based on the transfer of chiral edge states present 
in the populated sector (+) = {g^^e^} system, to the empty sector (— ) = {^|, e^}. A crucial point is to 
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Fig. 8. (Color online) Shelving method to detect edge states. The scheme is presented for ^^^Yb, a two-electron 
atom, with g the electronic ground state manifold and e a metastable excited state manifold, both with spin 
1/2. (a) The laser G couples the states and e^, but also g^ and e^, to generate the artificial magnetic flux in 
both (+) and ( — ) sectors (cf. main text). The states g^ and are initially populated, while the other states are 
empty, (b) In order to probe the (+) = {g^,e^} system, one introduces a weaker additional laser P, coupling 
g^ ^ e^. The circles indicate relative populations in each state after the probe pulse, (c) After the probe pulse, 
atoms in the g manifold are dispatched using an auxiliary imaging transition g ^ f, where / designates a 
second excited state with short radiative lifetime, suitable for applying a strong radiation pressure pulse and 
kick the g atoms out of the trap, (d) The atoms are subsequently brought down to the g manifold using, 
e.g., adiabatic-passage techniques to ensure near-unit transfer efficiency, (e) The atoms are finally detected, 
without stray contributions from unperturbed atoms in other internal states, (f) The local density that would 
be imaged through the Shelving method. Here, the excitation fraction is N[q,ujL) — ^ p-{m,n). 



ensure that the topological edge states have the same structure in the initial and final states, so that a 
specific chirality is probed (cf. Section 4.2). To this end, the initially unpopulated states (— ) are also 
coupled by a laser generating the same gauge field as for (+). The degeneracies are split by a relatively 
strong magnetic field [48], AEa- = —gaTUiB^ where a = e^g denotes the ground or excited manifold, 
i ="[ / I, rrii = ±1/2 the nuclear spin quantum number, Qg/h ^ —750 Hz/G, and ge/h « —1250 Hz/G. 
A bias field B ~ 100 G thus leads to Zeeman shifts ~ ±25 kHz on the tt transitions and ^ ±100 kHz on 
the cr^ transitions: these shifts are large compared to typical Rabi frequencies of both the gauge-field 
and probe lasers (~ 1 kHz or less), and the effect of the different lasers can thus be treated independently. 

In order to probe the (±) system, one introduces a weaker additional laser, coupling g^ ^ e^. Due 
to the gauge coupling, a population will build up in the g^ state as well (roughly equal to that in the 
state since i? <C J/h). Those atoms will be missing in the final detection step. After probing, the lattice 
sites are isolated by rapidly raising the lattice height and switching off the artificial gauge field. Atoms in 
the g manifold are dispatched (possibly detected) using an auxiliary strong transition g ^ f. A natural 
choice for ^^^Yb is /=^Pi, with a linewidth 7//27r ~ 28 MHz much larger than any Zeeman splitting in 
g or e (thus prohibiting independent detection of atoms depending on their spin i / t)- Crucially, atoms 
in the e manifold are not in resonance with the imaging light and are therefore unaffected. The atoms 
are subsequently brought down to the g manifold using, e.g.^ adiabatic passage techniques, leaving the 
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state unaffected. A further imaging pulse allows to detect those atoms, initially excited by the probe 
pulse. One might worry that a fraction of atoms from the state could end up being transferred too, 
thus contaminating the final edge signal. Fortunately, the off-resonant excitation rate to "wrong" states 
will be smaller than the resonant rate by a factor scaling as ~ (i7/Z\^)^, with Az ^ \ge — 9g\B/2 a 
typical Zeeman splitting. Taking for example the parameters given in [43], one has J/h ^ 100 Hz, and 
f2 ~ 0.05J/h ~ 27r X 20 Hz, making the final contamination of by negligible (~ 10~^). 



(a) (b) (c) 




Fig. 9. (a)-(b) Density p(^_^{m,n) for q — ±4, 7 = 10 and ujl — 0.14J//i. The probe and system parameters in 
(a)-(b) are the same as in Fig. 7. (c) Excitation fraction versus probe frequency, for several potential shapes. 



4.2 Analysis of the Shelving method 

In order to study the effect of the Shelving method, we consider a simplified level scheme with two internal 
states only, denoted by the indices (±). We suppose that only the (+) sector is initially populated. The 
spatial profile of the coupling laser /l(^) is similar to the one used for Bragg excitations, but now the 
Pauli principle does not restrict the available final states, since the state (— ) is initially unoccupied. We 
write the coupling to the probe as 

where the operator cj^(^) creates a particle of the (±) sector in the eigenstate IV^^), and where /^^ has 

the same definition as in Eq. (9), since = = fiQ. The (+) sector is initially populated, such 
that the initial and excited states have the following forms 



|0) = |l...l I 0...0)(+)|0...0)(_), (19) 




where we suppose that f2 <C J/h to neglect higher order excitations. We note that k is no longer 
restricted by the Pauli principle, such that Uki = {ck — ei)/h may now take negative values. We follow 
the same treatment as for the Bragg scheme and we obtain the excitation fraction as 

N{q,LOL) = 2nf2H ^ ^ \IlfS^'\uu - ujl), (21) 

l<Ep k 
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which differs from Eq. (16) by the fact that the final states k are now unrestricted. However, we stress 
that the sum over the initial states, X1k£;f (^-'-)' ^^^^^ restricted by the Pauli principle: for 

00 L ^ J/^ and when Ep = —1.5 J, this allows to probe the edge states that are located in the first 
bulk gap only. This important fact leads to the asymmetry highlighted in the corresponding excited 
fraction N{q^uoL)^ illustrated in Fig. 9(a)-(b), which demonstrates the specific chirality of the extracted 
edge states. The excitation fraction N^q^uj^) — N^—q^uj^) is represented in Fig. 4(c), showing a clear 
resonance peak at low frequencies Hul <^ J. Interestingly, this result shows that the low-energy regime 
is still governed by the chiral edge states located in the bulk gap, although transitions are now allowed 
for all the states below E-p^ including the bulk states. Indeed, the signal N{q^ujL) — N{—q^ujL) remains 
small and flat in the "edge-bulk" region, while the chiral "edge-edge" peak stands even clearer than 
in the Bragg case (since more "edge-edge" transitions are allowed between states of same chirality). 
By setting the probe parameters {q^uoL) close to a resonance peak, one can now populate edge states 
into the (— ) sector and directly visualize them using state- selective imaging. The corresponding density 
/?(_)(m,n) = (lZ^(t)|n(_)(m, n)|lZ^(t)) is illustrated in Figs. 9(a)-(b) for q = ±4. The clear difference be- 
tween the two images, obtained with different signs of q but otherwise identical setups, is a direct proof 
of the chiral nature of the edge excitations populated by the probe. 

We finally stress that the condition 

is necessary in order to probe the edge-state structure. Indeed, if we consider a simpler scheme in which 
the (— ) sector is no longer subjected to a synthetic gauge potential, we find that N{q^uoL) — N{—q^ujL) ^ 
0. This observation shows that our scheme requires that the edge states of the (— ) sector should have 
the same chirality than the initially populated edge-states of the (+) sector, i.e. both systems should 
be subjected to the same synthetic magnetic flux. 



5 Conclusion 

In this work, we showed that cold atoms trapped in optical lattices and subjected to synthetic magnetic 
fields offer a unique platform to investigate the physics of topological edge states. In particular, we 
showed that an elegant detection scheme, refered to as the Shelving method, allows to identify and 
extract topological edge states in a highly controllable way, but also to image them using available 
imaging technics. 

In this sense, a cold-atom simulation of the quantum Hall setup, together with our detection scheme, 
offers the possibility to directly visualize the edge states and to access their dispersion relation, without 
performing any transport measurements. In particular, such an experiment could be realized in the 
interacting regime, where the dispersion relation of fractional QH edge states could be accessed. 
Let us stress that our detection method does not rely on the lattice, nor on the setup which generates 
the synthetic magnetic field (e.g. laser-induced [39], lattice shaking [41], atom-chip [20], rotation, ...). 
Consequently, our method applies for a wide range of cold-atom realization of 2D topological phases, 
displaying a circular geometry. In particular, exploiting atomic species with many internal states, one 
could easily extend the present detection scheme to probe the physics of helical edge states exhibiting 
the quantum spin Hall effect [9,20,49,50]. 
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